arXiv: 1504.01787vl [cond-mat.quant-gas] 8 Apr 2015 


Collision-induced frequency shifts in bright matter-wave solitons and soliton molecules 


A. D. Martin 

Department of Physics, University of Otago, Dunedin 
(Dated: April 9, 2015) 

A recent experiment has detected collision-induced frequency shifts in bright matter-wave solitons for the 
first time [Nat. Phys. 10, 918 (2014)]. Using a particle model, we derive the frequency shift for two solitons in 
a harmonic trap, and compare it to the results of the recent experiment and reported theoretical approximation. 
We find regimes where the frequency shift is much smaller than previously predicted, and propose experiments 
to test these findings. We also predict that reducing the experimental trap frequency will reveal for the first time 
soliton molecules or soliton bound states in a cold atoms system. The dynamics of such bound states are found 
to be both highly phase-dependent, and sensitive to the residual 3D nature of the experiment. 

PACS numbers: 03.75.Lm, 


Since the first realisations of bright solitons in Bose- 
Einstein condensates (BECs) nua, many uses for solitons 
have been suggested including as surface-probes 13 m, for 
Bell-state generation HD and for interferometry The 

longstanding interest in bright solitons has recently been re¬ 
vived by an experiment on soliton collisions ifTOl . 

Solitons are solitary waves robust with respect to collisions; 
they survive these collisions unscathed up to a shift in their po¬ 
sitions and phases. However, the vital characteristic of bright 
solitons - the collisional shifts - has not been observed in BEC 
until recently. The experiment of Rice University nni of two 
solitons repeatedly colliding in a harmonic trap, demonstrated 
these shifts for the first time. The position shifts were revealed 
by their effect on the solitons’ frequency of oscillation in the 
trap. The authors of Ref. iTTOl compared their results with an 
approximate theoretical curve, which agreed reasonably well 
with the experiment in the tested regime. 

In this Letter we provide an improved prediction of the fre¬ 
quency shift based on a particle model that reproduces the ex¬ 
act position shifts in the untrapped case lfTTHT3]| . We identify 
regimes in which there is a measurable difference between our 
predictions and the theory of Ref. ifTOl and propose simple 
modifications to the experiment to explore these regimes. We 
also suggest modifications which will produce exotic soliton 
molecules (bound states) 0 Optical soliton molecules have 
been produced experimentally lfT4l . but until now such states 
have not been produced in a cold atoms system. We predict 
that a reduction in the trap frequency of the experiment ITOl . 
along with modifications to control the relative phase, would 
produce soliton molecules in BEC for the first time. We in¬ 
vestigate the 3-dimensional (3D) corrections to the dynamics, 
which affect both the stability and the mean-field dynamics of 
the BEC, particularly for the soliton molecules. 

Bose-Einstein condensates with tight radial trapping are in 
most circumstances well-described by the 1-dimensional (ID) 
Gross-Pitaevskii Equation (GPE) |[T5]| : 

d mo/x^ 2 

= ^— i// + giDN\ifrfi//, ( 1 ) 

where giD = IfiojrCis^ ^ is the total atom number, oj and ojr 
are the axial and radial trapping angular frequencies, and m 


and Os are the mass and s-wave scattering length of the atomic 
species. In experiments, Us may be tuned to vary the sign and 
magnitude of the interactions, e.g., negative Og permits bright 
soliton solutions ifTbl . We note that the ID GPE breaks down 
in regimes where quantum or thermal fluctuations ifTTll and/or 
3D effects become significant flEli . Quantum and thermal ef¬ 
fects are not expected to play a large role in the situations 
considered in this paper. However, the slightly 3D nature of 
the geometry imposes a limit on the interaction strength, be¬ 
yond which the condensate will collapse lfT0l[T9l . This limit 
is characterised by a critical atom number \Nc\. The param¬ 
eter Nc = 0.61 Or I as, where Oy is the radial harmonic length 
of the trapping potential, is also used here and in Ref. cni 
to characterise the strength of the interactions, since Ng/Nc is 
proportional to gwN, and Ng ^ N/2 is the number of par- 
ticle^er soliton. Adding the quintic nonlinearity g 2 W^^ to 
Eq. ([T), where g 2 = 24In (413)N^hcOrO^^, yields a ID equation 
which approximates the integrability-broken nature of the full 
3D GPE 1201 . We use the GPE equation modified with this 
quintic term to verify the regimes of validity of our results 
based on Eq. Q- 

In previous work |[TT1 - IT^ . the following soliton-soliton in¬ 
teraction potential was employed to model solitons as individ¬ 
ual particles: 


y = -2tjitj2 (tji + 772 ) sech ^ 


2mm ~ 
—I—^ 

m + ^2 


( 2 ) 


where q = m {gwl Nqlfi^ is the dimensionless form of the rel¬ 
ative position coordinate, q, of the solitons, and qt are effec¬ 
tive soliton mass parameters, which take the value q = 1/8 
for two solitons of equal size (121 [13. This potential pro¬ 
duces the exact position shifts for solitons emerging from col¬ 
lisions within Eq. Q when the axial potential oj = 0. The 
interaction potential was shown numerically also to provide 
a good description when combined with an external trapping 
potential in regimes where the solitons were well-separated 
between collisions, particularly when the solitons collided in- 
phase (13 . 

Within the particle model, the oscillation frequency of the 
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solitons’ relative position dynamics is given by 


^ = 2711 


X qo 

qo 


dx 

^{oj/4) {ql - x^) - Arf (sech^ (;/§o) - sech^ {rjx)) 


(3) 

Here, +§o are the turning points of the dimensionless relative 
position coordinate, q. For solitons initially at rest (as in the 
recent experiment Ool and all cases analysed in this paper), 
is equal to the initial soliton separation. The frequency-shift 
may then be written Au; = Q.- oj. We have found analytical 
solutions to Eq. ([^, in cases where the solitons are initially 
close or far apart. In the limit q{)m\giD\Njfi 1, i.e., where 

the solitons are at all times extremely close, we find: 



^(JJ + l6\giDf 

OJ OJ 


(4) 


For large initial separation, where qomlgioW/% » 1 and 
ff ln(l + l6Tf-\gw\^N^/lfo/ql) /27Tqoqm\gw\N 1, we can 
assume that the position-shifts generated by Eq. © occur on a 
timescale much faster than the trap period, and the oscillation 
period is found by merely adding the time-shifts generated by 
Eq. Q during two collisions to the trap period, yielding: 
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The behaviour of the frequency shift between these two 
limits can be explained qualitatively by considering the effec¬ 
tive potential, illustrated in Fig. The potential has a clear 
bi-modal character, with a narrower interaction mode within 
the harmonic potential. In the limit where q{)m\gij)\Njfi 1, 
the solitons are extremely close with respect to the width 
of the inter-soliton potential (which becomes large when the 
solitons’ interaction strength decreases). The dynamics are 
completely described by small harmonic oscillations, which 
increase in frequency as the interaction strength increases. 
When the interactions are increased further, or the solitons’ 
initial separations are increased, the solitons may still be 
strongly bound within the interaction potential, but this po¬ 
tential is no longer effectively harmonic, and the frequency in¬ 
crease with interaction strength is less rapid. For even stronger 
interactions or wider separations, the inter-soliton potential 
becomes narrow compared with the initial soliton separation, 
and the solitons effectively ‘escape’ the inter-soliton potential 
each collision. In this regime, the frequency shift levels off, 
and starts to decrease as the width of the inter-soliton poten¬ 
tial vanishes in comparison with q^. 

We provide for comparison the theoretical approximation 
derived in Ref. Co) based on the interaction Hamiltonian of a 
linear superposition of two Gaussian wavepackets: 


^(JJ _ -gipNal 
27iq^Q%(jij 


(7) 



FIG. 1. (Colour online.) Effective potential for the system of two 
solitons in a harmonic trap, with scattering length = -0.57ao and 
N = 56000, and the mass parameter corresponding to the Lithium- 
7 atomic species. The dashed curved shows the harmonic trapping 
potential alone. The trap frequencies are (a) io -27: x 31Hz, and (b) 
(jj - 2?: X 3Hz. The solid (red) vertical lines are at q^ - 26yL/m, and 
the dashed (red) vertical lines at q^ - 106yum. 


where is the axial harmonic length of the trapping poten¬ 
tial. As opposed to the behaviour of the dynamics in the par¬ 
ticle model, within Eq. ([T]), the shift is linear in the interaction 
strength. The most marked difference between the prediction 
of the particle model and that of Eq. 0 is that the frequency 
shift in Eq. (|7]) diverges to infinity, rather than to zero, as the 
interaction strength goes to infinity. In the opposite limit the 
particle model goes to zero quadratically, rather than linearly. 
However, the important comparison is for experimentally ac¬ 
cessible regions (NsINc > -1), and we are most interested 
in regimes where the difference between the predictions is 
large enough to resolve in experiment. In Ref. ifTOl . the initial 
soliton separation, q^ = 26/im; number of atoms per soliton. 
Ns = N/2 = 28000; and radial and axial trap angular frequen¬ 
cies, (jjy - 271 X 254Hz and oj = 27 i x 31Hz. The scattering 
length Us was varied between runs such that NsINc took values 
between -0.53 and -fO.55. Note that in the regimes of posi¬ 
tive Ns/Nc, the wavepackets are not strictly solitons, and such 
regimes are not considered in this paper. 

Figure shows the frequency shift from numerical evalu- 
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FIG. 2. (Colour online.) Relative frequency shifts versus NJNc for (a), (c) (jj - In Y. 31Hz, and (b), (d) cjl> = 2n x 3Hz. Plots (a), (b) are for 
initial soliton separation qo = 26yL/m, and (c), (d) are for qo - 106yL/m. The solid lines were determined by numerical evaluations of Eq. 

The dot-dashed lines are the asymptotic limit of small separations and weak interactions, given by Eq. 0- the dotted lines are the asymptotic 
limit of large separations and strong interactions given by Eq. The dashed (red) lines give the approximate result from Eq. 0 and (blue) 
bars give the experimental results Co). 


ation of Eq. along with the curves for Eqs. 0 and 0 
representing the limiting behaviour, and the theoretical ap¬ 
proximation provided in Ref. ifTOl . The relevant region of 
the experimental parameters lies within Fig.j^a); while Figs. 
[^b)-(d) show some other regimes of interest - with different 
initial soliton separations or trap frequencies than in the ex¬ 
periment. The main deviation between the particle model and 
Eq. ([^ is obvious in the large negative Ns/Nc limit, where the 
frequency in the particle model levels off and then starts to de¬ 
crease. This discrepancy is particularly apparent in Figs.j^b) 
and (d), where the axial trapping frequency is a factor of ap¬ 
proximately 10 less than in the recent experiment. However, 
for the parameters of recent experiment [Fig.j^a)] the differ¬ 
ence between the curves is probably within experimental er¬ 
ror; for large initial soliton separation and tight trapping [Fig. 
[^c)], the agreement between the curves is extremely close. 

Before issuing confident experimental predictions from the 
above particle model results, there are two reasons to verify 
the model’s performance. Firstly, the particle model (and also 


the theoretical approximation of Ref. HOl ) assumes that the 
solitons separate between collisions, and approach at a suffi¬ 
cient speed such that the relative-phase between the solitons 
plays no effect. It is clear from the results of Ref. |[T3]| that 
for solitons slow enough that the collision time approaches 
the order of the trap period, out-of-phase collisions will have 
smaller frequency-shifts than those predicted by the particle 
model, whereas in-phase collisions are still well described. In 
the regimes where the initial separation was qo = 26jim cni, 
the ansatz used to derive the particle model consists of initially 
slightly overlapping solitons, even for the narrowest solitons 
considered in the experiment. This obviously contradicts the 
assumption that the solitons always separate completely be¬ 
tween collisions. Secondly, in some regimes, 3D effects to 
lead to condensate collapse during collisions, especially for 
^s/^c > “0-5 and in-phase solitons ifTOlfT^ . 

We first evaluate on the phase dependence of the position 
shifts by integrating the GPE without the quintic term. Fig. 
I^a)-(b) shows GPE simulations two solitons for the experi- 
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FIG. 3. (Colour online.) Atomic density within the GPE for initial 
separation: (a)-(d) qo - 26jum, (e)-(f) qo = axial trap angular 

frequency: (a)-(b) cx) - In x 31Hz, (c)-(f) to - 2n x 3Hz. The 
relative phase between the solitons is (a), (c), (e) 0 = 0, (b), (d), (f) 
(f) = 71. The scattering length = -0.51 uq, the mass is that of the 
Lithium-7 atomic species and N - 56000 such that NJNc = -0.53. 
Trajectories within the particle model are given by full (red) lines, 
and the trajectories of non-interacting particles are given by dashed 
(red) lines for comparison. 


mental parameters with Ns/Nc = -0.53 for in- and out-of- 
phase collisions. The in-phase collisions are well described by 
the particle model, and, as expected, the predicted frequency 
shift is too large for the out-of-phase collisions. Interestingly, 
we find that when the axial trap is weakened to oj = 27rx 3Hz 
[Fig.[^c)-(d)], the in-phase solitons form a bound state (soli- 
ton molecule) which is surprisingly well-described by the par¬ 
ticle model (until it eventually de-phases at longer timescales). 
The out-of-phase solitons (which seem to repel each other) are 
badly described by the particle model. Similar solitons with 
relative phases other than 0 and tt will have intermediate be¬ 
haviour. We expect that the theoretical model employed in 
Ref. [171 would provide a better description of the GPE dy¬ 
namics of these soliton molecules if it could be extended to 
incorporate the (weak) external trapping potential. 

For larger initial soliton separations than used in the recent 
experiment, the particle model agrees very well with the GPE 
for any value of the relative phase. We illustrate this by sim¬ 
ulating the system for an initial separation of four times that 
of recent experiment: qo = 106yum. For the experimental trap 
frequency, the frequency shift is too small to detect practically. 
We instead consider a weaker trap frequency of co = 271 x 3Hz, 
resulting in slower solitons with easily observable frequency 
shifts. The agreement between the model and the GPE is il¬ 




FIG. 4. (Colour online.) (a)-(b) Atomic density within the GPE for a 
bound state without (a) and with (b) a quintic nonlinearity simulating 
the 3D nature of the EEC. Here, NJNc = -0.41, oj = 21 pi x 3Hz. 
(c)-(d) Atomic density within the GPE for a bound state with a quin¬ 
tic nonlinearity for in-phase (c) and out-of-phase (d) collisions. Tra¬ 
jectories within the particle model are given by full (red) lines, and 
the trajectories of non-interacting particles are given by dashed (red) 
lines. 


lustrated in Fig. [^e)-(f). Note that it is this regime in which 
there is greatest improvement of the particle model over Eq. 

0 . 

We now investigate 3D effects by integrating the ID GPE 
with the quintic nonlinearity. We find that for interaction 
strength Ng/Nc = -0.53, the in-phase collisions are unstable, 
to collapse, but not the out-of-phase collisions [see Fig. [^c)- 
(d)]. This suggests that the in-phase collisions of Fig.[^a),(c) 
and (e) would not be realisable. However, slightly reducing 
the interaction strength to Ng/Nc = -0.41 produced qualita¬ 
tively similar states which do not collapse. We find that the 
quintic term increases the oscillation frequency by a factor of 
approximately 1.5 [see figure 4(a)-(b)]. The frequency shifts 
of the unbound states remain unaffected [see Fig.j^d)], and it 
is here that the particle model is most useful. 

In conclusion, we demonstrated the existence of regimes 
described better by the particle model than the theoretical ap¬ 
proximation Cni, i.e., those of weak axial trap frequency and 
stronger interactions. We propose the extension of the recent 
experiment cni to verify the particle model predictions. We 
also suggest that soliton molecules can be created if greater 
control of the soliton phase is possible, e.g., by applying a 
light-sheet potential to half of the condensate ll2TH23l . How¬ 
ever, such states are highly dependent on relative-phase, and 
also on 3D effects. In future work we propose extending the 
formalism of Q to include the effects of a harmonic trapping 
potential and 3D effects in order to describe soliton molecules. 

We thank D. Schumayer and R B. Blakie for useful discus¬ 
sions, and P. B. Blakie for assistance with preparation of Fig, 
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